Plot_Sel <- function (dat=msm1,ylim=c(0,2.5)) {
  par(mfrow=c(1,3))
  plot(dat$srv_sel_1,ylim=ylim,main="Pollock",ylab="Selectivity")
  plot(dat$srv_sel_2,ylim=ylim,main="Pcod",ylab="")
  plot(dat$srv_sel_3,ylim=ylim,main="ATF",ylab="")
  plot(dat$fsh_sel_1,ylim=ylim,main="Pollock",ylab="Selectivity")
  plot(dat$fsh_sel_2,ylim=ylim,main="Pcod",ylab="")
  plot(dat$fsh_sel_3,ylim=ylim,main="ATF",ylab="")
}



